部分均匀荷载作用下圆板塑性极限统一解

本文采用俞茂宏双剪统一屈服准则求解部分均匀荷载下简支圆板的塑性极限，已有的最大主应力准则、Ｔｒｅｓｃａ准则，Ｍｉｓｅｓ准则的解答均是该解答的特例或线性逼近。该解答可以适合于多种工程材料，具有普遍性和适用性，从得出的解析解和图示中可得出一些新的有益的认识     

沿垂等温平板的自然对流问题的近似解析解

本文用加权残值得到了沿垂等温平板的自然对流问题的近似解析解。和已有准确解的比较表明，在平板表面附近，近似解有很好的精度，而且速度场的精度比温度场更好一些。文末还对配域法的准收敛问题作了说明。     

用LDV系统测量转盘萃取塔内流场

本文应用激光多普勒测速仪对二种转盘萃取塔模型进行切向速度、轴向速度和湍流度测量。实验结果表明，在ＲＤＣ是切向速度比较大，因而有较大的搅拌效果，它的内部流场比较均匀稳定，有利于传质，所以ＲＤＣ是广泛应用于化学工业中的一种液－液萃取设备。     

Singularity-free augmented Lagrangian algorithms for constrained multibody dynamics

After a general review of the methods currently available for the dynamics of constrained multibody systems in the context of numerical efficiency and ability to solve the differential equations of motion in singular positions, we examine the acceleration based augmented Lagrangian formulations, and propose a new one for holonomic and non-holonomic systems that is based on the canonical equations of Hamilton. This new one proves to be more stable and accurate that the acceleration based counterpart under repetitive singular positions. The proposed algorithms are numerically efficient, can use standard conditionally stable numerical integrators and do not fail in singular positions, as the classical formulations do. The reason for the numerical efficiency and better behavior under singularities relies on the fact that the leading matrix of the resultant system of ODEs is sparse, symmetric, positive definite, and its rank is independent of that of the Jacobian of the constraint equations. The latter fact makes the proposed method particularly suitable for singular configurations.     

Resolving Controversies in the Application of the Method of Multiple Scales and the Generalized Method of Averaging

I compare application of the method of multiple scales with reconstitution and the generalized method of averaging for determining higher-order approximations of three single-degree-of-freedom systems and a two-degree-of-freedom system. Three implementations of the method of multiple scales are considered, namely, application of the method to the system equations expressed as second-order equations, as first-order equations, and in complex-variable form. I show that all of these methods produce the same modulation equations.I address the problem of determining higher-order approximate solutions of the Duffing equation in the case of primary resonance. I show that the conclusions of Rahman and Burton that the method of multiple scales, the generalized method of averaging, and Lie series and transforms might lead to incorrect results, in that spurious solutions occur and the obtained frequency–response curves bear little resemblance to the actual response, is the result of their using parameter values for which the neglected terms are the same order as the retained terms. I show also that spurious solutions cannot be avoided, in general, in any consistent expansion and their presence does not constitute a limitation of the methods. In particular, I show that, for the Duffing equation, the second-order frequency–response equation does not possess spurious solutions for the case of hardening nonlinearity, but possesses spurious solutions for the case of softening nonlinearity. For sufficiently small nonlinearity, the spurious solutions are far removed from the actual response. But as the strength of the nonlinearity increases, these solutions move closer to the backbone and eventually distort it. This is not a drawback of the perturbation methods but an indication of an application of the analysis for parameter values outside the range of validity of the expansion.Also, I address the problem of obtaining non-Hamiltonian modulation equations in the application of the method of multiple scales to multi-degree-of-freedom Hamiltonian systems written as second-order equations in time and how this problem can be overcome by attacking the state-space form of the governing equations. Moreover, I show that application of a variation of the method of Rahman and Burton to multi-degree-of-freedom systems leads to results that do not agree with those obtained with the generalized method of averaging.Contributed by Prof. R.A. Ibrahim.     

汽轮机转子在气流力和油膜力作用下的非线性动力学特性

为了分析转子在油膜力和气流激振力共同作用下的非线性振动特性，本文以短轴承支撑的不平衡刚性对称Jeffcott转子系统为研究对象，首先分析转子在非稳态油膜力作用下的振动特性，然后分析转子在油膜力和气流激振力共同作用下的非线性振动特性。采用数值模拟的方法研究了系统的分岔和混沌特性，计算结果表明，考虑气流激振力和油膜力共同作用下的转子系统与仅考虑油膜力的转子系统相比，在相对进气速度v=30m／s时，随着无量纲转速ω的增加。二者都出现了周期运动和混沌运动多次交替出现的复杂运动特性，但是前者首次出现倍周期分岔和混沌运动时的转运提前，在定转速情况下，随着v的增大，系统最终在经历周期运动之后进入混沌运动，而且圆盘中心的最大振幅随着v的增大而增大。     

无网格算法在多段翼型流动计算中的应用

研究了一种求解欧拉方程的无网格算法，发展出了一套布点及点云自动生成的方法；在点云离散的基础上，采用最小二乘法求解矛盾方程的方法来求取空间导数，进而获得数值通量；采用四步龙格-库塔方法进行时间推进，并引入当地时间步长和残值光顺等加速收敛措施。通过对NA-CA0012翼型的跨音速流动和多段翼型复杂绕流的数值模拟，验证了上述无网格算法的正确性和实用性。     

Deformation model for brittle materials and the structure of failure waves

Constitutive equations that describe the experimentally observed failure waves are proposed to model inelastic strains of brittle materials. The complete system of equations is hyperbolic, each equation of this system has divergent form. The model is based on the assumption that continual failure is the process of transition from an intact state to a “fully damaged” state described by the kinetics of the order parameter. The structure of stationary traveling compressive waves is analyzed using a simplified model. It is shown that in a certain range of amplitudes, the wave splits into an elastic precursor and a failure wave. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 3, pp. 164–172, May–June, 2007.     

On the basic laws of anisotropic viscoelasticity

The main aim of this work is to develop a consistent formulation of the rheological behavior for different anisotropic polymer systems. The unified theory of anisotropic viscoelasticity is developed based on the symmetry principles. The Maxwell rheological equation is extended to nonsymmetric anisotropic liquids. Transitions from the most general anisotropy to particular cases of anisotropy are established. It appears that the coupled relaxation of symmetric and antisymmetric stresses is a natural phenomenon in nonsymmetric viscoelasticity. Within the concept of an internal state variable, a stress–order relation is derived for a fully nonlinear case. The order tensor dynamics is also considered. A simple method of deriving the equation of the internal rotational motion is developed for the general macroscopic anisotropy. This paper was presented at the 3rd Annual Rheology Conference, AERC 2006, April 27–29, 2006, Crete, Greece     

方形管内楔形障碍物对火焰结构与传播的影响

通过实验与数值模拟方法对CH4/空气预混火焰在有楔形障碍物的卧式燃烧方管内的传播进行了研究。采用多镜头Cranz Schardin高速摄像机和压力传感器等实验设备获得了高清晰度的障碍物诱导火焰失稳的分幅时序照片以及障碍物背风表面压力变化曲线。数值模拟则基于RANS方法与EDU-Arrhenius燃烧模型,计算结果与实验结果基本相符,反映了火焰在管内传播与变形的详细过程。通过综合分析实验与计算结果,得到了由楔形障碍物导致的火焰加速与变形的内在机理,揭示了火焰传播过程中由层流燃烧向湍流燃烧转捩的本质。